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The possible exotic nuclear properties in the neutron-rich Ca, Ni, Zr, and Sn isotopes are examined 
with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of Green’s function 
method. The pairing correlation, the couplings with continuum, and the blocking effects for the 
unpaired nucleon in odd-A nuclei are properly treated. The Skyrme interaction SLy4 is adopted for 
the ph channel and the pairing interaction DDDI is adopted for the pp channel, which well reproduce 
the experimental two-neutron separation energies S2, and one-neutron separation energies Sy. It is 
found that the criteria S, > 0 predicts a neutron drip line with much smaller neutron numbers than 
that of Son > 0. Due to the the unpaired odd neutron, neutron pairing energies —Epair in odd-A 
nuclei are much smaller than the neighboring even-even nuclei. By investigating the single-particle 
structures, the possible halo structures in the neutron-rich Ca, Ni, and Sn isotopes are predicted, 
where a sharp increase of rms radii with significant deviations from the traditional r œ A‘? rule 
and diffuse spatial distributions in densities are observed. By analyzing the contributions of various 
partial waves to the total neutron density pz; (r)/p(r), the orbitals located around the Fermi surface 
especially those with low angular momenta play significant roles in the extended nuclear density 
and large rms radii. The number of the neutrons N) (No) occupying above the Fermi surface 
An (continuum threshold) are discussed, the evolution of which as a function of the mass number A 
in each isotopes consist with those of pairing energy, supporting the key role of pairing correlation 
in halo phenomena. 


PACS numbers: 21.60.-n, 21.10.Gv, 21.10.-k, 21.60.Jz 


I. INTRODUCTION neutron or proton Fermi surfaces are usually very close to 
the continuum threshold. With the effects of the pairing 


correlation, the valence nucleons have a certain probabil- 


The study of exotic nuclei far from the 6 stability line 
is a very challenging frontier topic in nuclear physics ex- 
perimentally [1—4] and theoretically [5-10]. The unstable 
nuclei with extreme N/Z ratios, which are weakly bound 
systems, have shown many exotic properties being differ- 
ent from the stable nuclei, such as the halo structures [11— 
17], the changes of traditional magic numbers [18-23], 
and new nuclear excitation modes [24, 25], which may 
also herald new physics. The study of exotic nuclei is 
not only crucial for the fully understanding of the rich 
nuclear structures and properties, but also important for 
the investigations of element synthesis and nuclear astro- 
physics [26]. However, due to the extreme short lifespan, 
the cross sections for synthesizing exotic nuclei are ex- 
tremely small, which makes it extremely difficult to cre- 
ate them experimentally. As a result, more and more 
advanced large-scale radioactive beam facilities and up- 
dated detector techniques have been built, upgraded, or 
planned in the worldwide [27-34]. Meanwhile, abundant 
theoretical studies of exotic nuclei provide a useful way 
for the layout of experiments and analysis of experimen- 
tal results [35-37]. 


In the exotic nuclei, especially the drip line nuclei, the 
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ity to be scattered into the continuum and occupy the 
resonant states therein, making the nuclear density dis- 
tributions very diffuse and extended. It is therefore es- 
sential to treat the pairing correlations and the couplings 
to continuum properly in the theoretical descriptions of 
exotic nuclei [38-42]. Besides, in the one-neutron halo 
nuclei, such as 3tNe [43] and 3’Mg [44], the blocking ef- 
fect [45] should also be considered to treat the unpaired 
odd nucleon. 


The Hartree-Fock-Bogliubov (HFB) theory has 
achieved great successes in describing exotic nuclei 
with a unified description of the mean field and pairing 
correlation via Bogoliubov transformation [45]. Different 
models based on HFB theory have been taken to study 
the exotic nuclei, such as the Gogny-HFB theory [46], 
the Skyrme-HFB theory [47], the relativistic continuum 
Hartree-Bogoliubov (RCHB) theory [39], and the density 
dependent relativistic Hartree-Fock-Bogoliubov (RHF- 
B) theory [48]. To explore the halo phenomena in 
deformed nuclei, these models have been extended to the 
deformed framework, such as the deformed relativistic 
Hartree-Bogoliubov (DRHB) theory [16, 49, 50] and the 
coordinate-space Skyrme-HFB approach [51-54]. 


Traditionally, these H(F)B equations are often solved 
in configuration spaces, i.e., the basis expansion 
method [55]. However, the calculations are strongly re- 


lated to the space size and the shape of the expanded 
basis. For the harmonic oscillator basis, which has a tail 
in a shape of Gaussian function, although it is very effi- 
cient to describe the stable nuclei, significant difficulties 
were encountered when applying to exotic nuclei. Bases 
with proper shapes, such as the Woods-Saxon basis [56] 
and the transformed harmonic oscillator basis [57, 58] are 
often taken for the exotic nuclei. For example, to explore 
deformed halos [59, 60], the DRHB theory based on a 
Woods-Saxon basis [16, 49] has been developed. In con- 
trast with the basis expansion method, solving the HFB 
equation in the coordinate space is believed to be more 
effective. In the coordinate space, discretized method 
with the box boundary condition has been widely used, 
by which a series of discrete quasiparticle levels can be 
obtained easily. However, some flaws have been iden- 
tified for this method, such as the nonphysical drop of 
the nuclear densities at the box boundary, the continua 
and resonant states are discretized, and some unphysi- 
cal states are included. In contrast, the Green’s func- 
tion (GF) method [61-63] in the coordinate space can 
avoid these problems and owns great advantages, which 
can describe the asymptotic behaviors of wave functions 
properly, provide the energies and widths of the resonant 
states directly, and treat the bound states and the con- 
tinua on the same footing. 


Due to these advantages, the Green’s function method 
has been applied extensively in nuclear physics to study 
the contribution of the continuum to the nuclear struc- 
tures and excitations. As early as in 1987, Belyaev et al. 
have constructed the Green’s function for the HFB equa- 
tion [64]. Afterwards, this HFB Green’s function was 
applied to the quasiparticle random-phase approxima- 
tion (QRPA) [65], which was further used to describe the 
collective excitations coupled to the continuum [66-71]. 
In 2009, the continuum HFB theory in a coupled channel 
representation was developed to explore the effects of the 
continuum and pairing correlation in deformed neutron- 
rich Mg isotopes [72]. In 2011, Zhang et al. developed the 
fully self-consistent continuum Skyrme-HFB theory with 
Green’s function method [73], which was further applied 
to investigate the giant halos [74] and the effects of pair- 
ing correlation on the quasiparticle resonances [75, 76]. 
In 2019, to explore the halo phenomena in the neutron- 
rich odd-A nuclei, the self-consistent continuum Skyrme- 
HFB theory was further extended by including the block- 
ing effect [77]. In recent years, Green’s function method 
has also been adopted for the covariant density function- 
al theory (CDFT) [78-86] in the studies of nuclear struc- 
ture. For example, by introducing the Green’s function 
method to the relativistic mean filed theory (GF-RMF), 
the single-particle level structures including the bound 
states and resonant states and the pseudospin symme- 
tries therein are investigated for neutrons [87, 88], pro- 
tons [89], and A hyperons [90]. Besides, it has been con- 
firmed that exact values of the energies and widths could 
be obtained by searching for the poles of Green’s func- 
tion or the extremes of the density of states in disregard 


of the widths of resonant states [91, 92]. By combining 
the Green’s function method with the RCHB theory, the 
pairing correlation and continuum are well described in 
the giant halos of the Zr isotopes [93]. By extending the 
GF-RMF model to the coupled channel representation, 
the halo candidate nucleus ?’Mg reported experimental- 
ly is analyzed and confirmed to be a p-wave one-neutron 
halo according to the Nilsson levels [94]. In addition, the 
complex-scaled Green’s function method [95] has been 
established as a powerful tool for the exploration of res- 
onant states, which was further extended to the frame- 
work of relativistic mean field [96] and the deformed nu- 
clei [97]. Besides, the RMF-CMR-GF approach has been 
developed by combining the complex momentum repre- 
sentation method with the Green’s function method in 
the relativistic mean-field framework to study the ha- 
lo structures in the neutron-rich nuclei [98]. All those 
works have proved the great successes of Green’s func- 
tion method in the descriptions of the continuum. 


In this paper, the neutron-rich Ca, Ni, Zr, and Sn iso- 
topes are investigated systematically taking the contin- 
uum Skyrme-HFB theory formulated with the Green’s 
function method in the coordinate space, which has treat- 
ed the pairing correlations, the couplings with continu- 
um, and the blocking effect for the odd unpaired nucleon 
properly. The paper is organized as follows: In Sec. II, 
we introduce briefly the continuum Skyrme-HFB theory. 
Numerical details are presented in Sec. III. After giving 
the results and discussions in Sec. IV, conclusions are 
drawn in Sec. V. 


II. THEORETICAL FRAMEWORK 


In the Hartree-Fock-Bogoliubov (HFB) theory [45], the 
pair correlated nuclear system is described in terms of 
independent quasiparticles by the Bogoliubov transfor- 
mation. The HFB equation in the coordinate space [38] 
is, 
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where F; is the quasiparticle energy, ¢;(ro) is the quasi- 
particle wave function, and is the Fermi energy deter- 
mined by constraining the expectation value of the nu- 
cleon number. The HF Hamiltonian h(ro,r’o’) and the 
pair Hamiltonian h(ro,r’o’) are obtained by the varia- 
tion of the total energy functional with respect to the par- 
ticle density p(ro, r'o’) and the pair density p(ro,r’o’), 
respectively. The solutions of the HFB equations have 
two symmetric branches, one is positive (Æ; > 0) with 
the quasiparticle wave function ¢;(ro) while the other 
is negative (—E; < 0) with the conjugate wave function 
$;(ra). The wave functions ¢;(ro) and ¢;(rc) are in two 


components and they are related as, 
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with y(ra) = —2o0~(r,—c). Here and hereafter, we fol- 
low the notations in Ref. [65] for convenience. 

For an even-even nucleus, the ground state |®g) is a 
quasiparticle vacuum with the particle density p(ro, r'o’) 
and pairing density p(ro,r’o’) determined by 
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p(ro,r’o') = 


p(ro,r’a') = 


where cl, and Cro are the particle creation and annihi- 
lation operators, respectively. The densities can be u- 
nified as a generalized density matrix R(ro,r’o’), with 
p(ro,r’a’) and p(ro,r’o’) being the “11” and “22” el- 
ement, respectively. With the quasiparticle wave func- 
tions, R(ro,r’o’) can be written in a simple form, 


= Daal (ro)o!(r'o"). (4) 


R(ro,r’o') 


For an odd-A nucleus, the last odd nucleon is unpaired, 
for which the blocking effect should be considered. The 
nuclear ground state in this case is a one-quasiparticle 
state |®;), which can be constructed based on a HFB 
vacuum |®9) as 


|®1) = £$ |p), (5) 


where Bh is the quasiparticle creation operator and ip 
denotes the blocked quasiparticle level occupied by the 
odd nucleon. Accordingly, the particle density p(ro, r'o’) 
and the pairing density 6(ro,r’a’) are 


p(ro,r’a') = (dilet 0 Cro|®1), (6a) 
pro, r'o’) = (1 |era’ Cro|®1), (6b) 
and the generalized density matrix R(ro, r'o’) becomes 
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where two more terms are introduced compared with 
those for the even-even nuclei. 

In the conventional Skyrme-HFB theory, the HFB e- 
quation (1) in the coordinate space is often solved with 
the box boundary condition, and a series of the dis- 
cretized eigensolutions including the quasiparticle en- 
ergy E; and the corresponding wave functions ¢,;(rc) 
could be obtained. Then the generalized density ma- 
trix R(ro,r’a’) can be calculated by summing those dis- 
cretized quasiparticle states as shown in Eqs. (4) and (7). 
We call this method the box-discretized approach. How- 
ever, the applicability of the box boundary condition in 


the description of exotic nuclei especially those close to 
the drip line is not good. A large enough coordinate s- 
pace (or box size) should be taken to describe the very 
extended density distribution. 

Green’s function method can avoid these problems 
of box-discretized approach, which impose the correc- 
t asymptotic behaviors on the wave functions, especial- 
ly for the weakly-bound sates and the continuum. The 
Green’s function G(ro,r’o’; E) with an arbitrary quasi- 
particle energy E defined for the coordinate-space HFB 
equation obeys, 
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which is a 2 x 2 matrix. The generalized density matrix 
R(ro,r’o’) in Eq. (7) can be calculated by the integrals 
of the Green’s function performed on the complex quasi- 
particle energy plane as 


R(ro,r'o’) = 2 $ dEG(ro,r'o'; E) 
Cr<o 
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(9) 


where the contour path Cp<o encloses all the negative 
quasiparticle energies —E; < 0, Cp encloses only the 
pole of —F;,, and Cr encloses only the pole of Fj,. 

In the spherical case, the quasiparticle wave functions 
ilro) and ¢;(ro) are only dependent on the radial parts, 
and they can be expanded as 


bilo) = Ton (Y jalo), (10) 
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where Yo) is the spin spherical harmonic, and 
Y; miro) = —20Y m (f — ø). Similarly, the general- 
ized density matrix R(ro,r’o’) and the Green’s function 
G(ro,r'o'; E) can be expanded as 


R(ro,r'o') =) Yh (Fo) Rij (1,7 )Y a(o"), (11a) 
ljm 
E 
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where Rı;(r,r’) and Gj;(r,r’; Æ) are the radial parts of 


the generalized density matrix and Green’s function, re- 
spectively. 

As a result, the radial local generalized density ma- 
trix R(r) = R(r,r) can be expressed by the radial HFB 
Green’s function Gi; (r, r’; E) as 


R(r) = Ð Ry(r,r) 
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From the radial generalized matrix R(r), one can easi- 
ly obtain the radial local particle density p(r) and pair 
density p(r), which are the “11” and “12” components 
of R(r), respectively. In the same way, one can express 
other radial local densities needed in the functional of 
the Skyrme interaction, such as the kinetic-energy den- 
sity T(r), the spin-orbit density J(r), etc., in terms of 
the radial Green’s function. For the construction of the 
Green’s function, see Refs. [73, 77]. 
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III. NUMERICAL DETAILS 


For the Skyrme interaction in the ph channel, the SLy4 
parameter set [100] is adopted. For the pairing interac- 
tion in the pp channel, a density dependent 6 interaction 
(DDDI) is taken, 


The pair Hamiltonian h(ro,r’o’) is then reduced to a 
local pair potential [47], 
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The strength of the pairing force Vo = —458.4 MeV-fm?, 
the density pọ = 0.08 fm~%, and other parameters 7 = 
0.71, œa = 0.59, which are constrained by reproducing 
the experimental neutron pairing gaps for the Sn iso- 
topes [68, 101, 102]. With these parameters, DDDI can 
reproduce the scattering length a = —18.5 fm in the !S 
channel of the bare nuclear force in the low density lim- 
it [101]. The cut-off of the quasiparticle states are taken 
with maximal angular momentum jmax = 25/2 and the 
maximal quasiparticle energy Ecut = 60 MeV. 

The HFB equation is solved in the coordinate space 


with the space size Rpox = 20 fm and mesh size dr = 
0.1 fm. To calculate the densities with the Green’s 
function, the integrals of the Green’s functions are per- 
formed along a contour path Cgeo, which is chosen 
as a rectangle with the height y = 0.1 MeV and the 
length Ecut = 60 MeV to enclose all the quasiparticle 
states with negative energies. For the odd-A nuclei, two 
more contour paths Cc and C, , which only enclose the 
blocked quasiparticle states at energies Fi, and —£j,, 
are introduced due to the blocking effect of the odd un- 
paired nucleon. For these details, one can refer to Re- 
f. [77]. To do the contour integration, an energy step of 
AE = 0.01 MeV on the contour path is adopted. 


IV. RESULTS AND DISCUSSION 


In Fig. 1, the two-neutron separation energies 
Son(N,Z) = E(N — 2, Z) — E(N, Z) are plotted for the 
even-even and odd-even Ca, Ni, Zr, and Sn isotopes. Red 
circles are those calculated by the continuum Skyrme- 
HFB theory with the SLy4 parameter set, in comparison 
with the results by the discretized method denoted by the 
blue triangles and the available experimental data [99] de- 
noted by the black squares. The differences between Son 
by the Green’s function method and discretized method 
are found slight. Good agreements with the experimental 
data could be observed, indicating the reliability of the 
continuum Skyrme-HFB theory in the prediction of neu- 
tron drip line. The traditional shell closures, i.e., N = 28 
in Ca isotopes, N = 50 in Ni isotopes, N = 50,82 in Zr 
isotopes, and N = 82 in Sn isotopes, could also be identi- 
fied, where the Sj, drops sharply. For example, in the Sn 
chain, Sən drops from 13.25 MeV at !2?Sn to 4.94 MeV 
at 134Sn with the neutron number exceeding the magic 
number N = 82. In the Ca, Ni, and Zr chains, the two- 
neutron separation energies quickly reach zero at large 
mass range, resulting in relative short neutron drip lines, 
which are °’Ca, ®°Ni, and !?3Zr, respectively. On the 
contrary, in the Sn chain, $2, keeps less than 1.0 MeV 
in a wide mass region after the gap of N = 82 and final- 
ly becomes negative until A = 178, suggesting !77Sn is 
a neutron drip line nucleus. Those weakly bound nuclei 
are quite interesting due to possible appearance of neu- 
tron halos although it is experimentally difficult to reach. 
In addition, the exploration of the neutron drip line and 
the determination on the limit of nuclear landscape are 
significantly important in nuclear physics. However, var- 
ious theoretical studies show that the predicated neutron 
drip line is very model dependent [103]. Moreover, dif- 
ferent physical quantities or criteria will also predict very 
different neutron drip lines. 

To explore the neutron drip lines in the Ca, Ni, Zr, 
and Sn isotopes, in Fig. 2 the single-neutron separation 
energies S,(N,Z) = E(N — 1,2) — E(N, Z) are plot- 
ted. Results by the continuum Skyrme-HFB theory with 
the SLy4 parameter set are denoted by red circles, which 
are consistent well with the experimental data [99] de- 
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FIG. 1: (Color online) Two-neutron separation energies S2n in the Ca, Ni, Zr, and Sn isotopes as a function of the mass number 
A calculated by the continuum Skyrme HFB theory with the SLy4 parameter set (filled red circles), in comparison with the 
results by the discretized method (open blue triangles) and the experimental data (filled black squares) [99]. 
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FIG. 2: (Color online) Single-neutron separation energies 
Sn as a function of the mass number A in the (a) Ca, 
(b) Ni, (c) Zr, and (d) Sn isotopes calculated by the continu- 
um Skyrme-HFB theory with the SLy4 parameter set (filled 
red circles), in comparison with the results by the dis- 
cretized method (open blue triangles) and the experimental 
data (filled black squares) [99]. 


noted by the black squares. Strong odd-even staggering 
is observed in all isotopes. In general, the S, in the 
even-even nucleus is around 2 ~ 3 MeV larger than the 
neighbouring odd-A nuclei, which is attributed to the un- 
paired odd neutron with vanishing pairing energy. As a 
result, compared with those in Fig. 1, the neutron drip 
lines determined by the one-neutron separation energy 
are greatly shortened. In the Ca, Ni, Zr, and Sn iso- 
topes, the drip line nuclei are °°Ca, ®®Ni, 1??Zr, and 
48Sn, respectively, the positions of which are indicat- 
ed by the black arrows. Outside the neutron drip line 


determined by Shn, the bound even-even nuclei behave as 
the interesting Borromean systems. For example, based 
on the bound nucleus ©°Ca, ©°Ca+n is unbound while 
Catn-+n is bound. 
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FIG. 3: (Color online) Neutron pairing energy —Epair as a 
function of the mass number A in the (a) Ca, (b) Ni, (c) Zr, 
and (d) Sn isotopes calculated by the continuum Skyrme-HFB 
theory with the SLy4 parameter set. The filled and open 
squares represent the results of the even-even and odd-even 
nuclei, respectively. 


In Fig. 3, we plot the neutron pairing energy Ebair, 
which is expressed as, 
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The red solid symbols are for the even-even nuclei and 
open symbols for the odd-A nuclei. The neutron pairing 
energies FEpair for the odd-A nuclei are much smaller than 
those of the neighboring even-even nuclei, due to the ab- 
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FIG. 4: (Color online) Single-particle levels for neutrons 
around the Fermi surface in the canonical basis for the (a) Ca, 
(b) Ni, (c) Zr, and (d) Sn isotopes calculated by the contin- 
uum Skyrme-HFB theory. The short dashed lines represent 
the Fermi surface. 


sent contribution of pairing energy by the unpaired neu- 
tron. This also explains why the drip line determined by 
the single-neutron separation energy Sn is much shorter 
than that obtained by the two-neutron separation ener- 
gy Soy. In addition, obvious shell effects are shown in 
the pairing energy. For example, in the Sn isotopes, zero 
pairing energies at N = 82 and N = 126 while reaching 
the largest value at half shell N = 102 could be observed. 
As a result, the traditional shell closures, i.e., N = 28, 40 
in Ca isotopes, N = 40,50 in Ni isotopes, N = 50,82 in 
Zr isotopes, and N = 82,126 in Sn isotopes, could also 
be obviously observed, which are consistent with those 
obtained in Fig. 1. In addition, a sub-shell N = 32 is 
also observed in Ca isotopes. 

In the following, the possible neutron halos in the Ca, 
Ni, Zr, and Sn isotopes are examined, especially in the 
weakly bound nuclei close to the neutron drip line, in 
which the Fermi surfaces are very close to the continu- 
um threshold and the valence neutrons could be easily 
scattered to the continuum. 

In Fig. 4, the neutron canonical single-particle struc- 
ture as functions of the mass number are plotted for the 
(a) Ca, (b) Ni, (c) Zr, and (d) Sn isotopes. For more 
details on obtaining the canonical single-particle level- 
s, please refer to Refs. [104, 105]. The neutron Fermi 
energy A, as well as the canonical single-particle ener- 
gies € are shown. As the neutron number increasing, 
the Fermi energy An in each chain is raised up and fi- 
nally reaches the continuum threshold while all the HF 
single-particle levels decreases. The traditional shell clo- 


sures, i.e., N = 28,40 in Ca isotopes, N = 40,50 in Ni 
isotopes, N = 82 in Zr isotopes, and N = 82 in Sn iso- 
topes, could be observed, where the big gaps could be ob- 
served obviously. Different single-particle structures are 
revealed in Ca, Ni, Zr, and Sn chains, based on which 
halos may be formed. In the Ni chain, above the shell 
closure of N = 50, there are several weakly bound states 
and low-lying positive canonical states in the continuum 
with small angular momenta, which favor the formation 
of halos. For example, in 8°Ni, around the Fermi sur- 
face, there are two weakly bound states, 2d5/2 and 381/2, 
and one low-lying positive canonical state 2d3/2. The Sn 
chain is quite similar to the Ni chain but has more ad- 
vantages for halo formation, where weakly bound states 
and the low-lying positive states in the continuum ex- 
ist above the N = 82 shell closure. Besides, the Fermi 
surface A, approaches to zero gradually and the Sn iso- 
topes in a large mass region are weakly bound. In the Ca 
chain, above the shell closure of N = 40, there is mainly 
the 1g9/2 state, which evolves from a canonical positive 
state in the continuum (A < 62) to a weakly bound lev- 
el (A > 64). Although there is a possibility for valence 
neutrons occupying 1gg/2 orbital, the contributed densi- 
ty is very localized due to the big central barrier. In the 
neutron-rich Ca isotopes, the low-lying positive canon- 
ical states in the continuum 31/2 and 2ds5/2 also play 
important roles. The formation of halos is most unlikely 
for Zr chain where the shell closure of N = 82 is located 
around the threshold of continuum and it’s hard for the 
valence neutrons to overcome the large gap and occupy 
the continuum. 


In Fig. 5, the neutron root-mean-square (rms) radii 


4 
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are plotted for the Ca, Ni, Zr, and Sn isotopes, which are 
obtained by the Skyrme-HFB calculations both with the 
Green’s function method (filled circles and solid lines) 
and box-discretized method (open triangles and dashed 
lines) employing the SLy4 parameter set. The radii 
r = bọ A!/3 in the traditional liquid-drop model (black 
lines) are indicated as well, where the coefficient bp could 
be determined by the radii of deeply bound nuclei. In 
the Ca, Ni, Zr, and Sn chains, they are r = 0.991A!/8, 
0.984A'/3, 0.957A!/3, and 0.961A!/%, respectively. In 
each chain, by adding more neutrons, the nuclear rm- 
s radii ryms increase steeply and deviate from the radii 
A‘/® rule. For example, in the Ca chain, compared with 
the isotopic trend in N < 20 with r = 0.991A!/3, the 
neutron rms radii in °°Ca and the heavier isotopes dis- 
play steep increases with N. In those mass area, pos- 
sible neutron halos may occur. Besides, obvious odd- 
even staggering of rms radii could be observed in the Sn 
chain, where the odd-A nuclei 4°!~!©°Sn own larger rms 
radii than the neighbouring two even-even nuclei. The 
detailed explanations could be found in Ref. [77]. The 
odd-even staggering phenomena in nuclear radii and nu- 
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FIG. 5: (Color online) Neutron root-mean-square (rms) radii 
Trms for the (a) Ca, (b) Ni, (c) Zr, and (d) Sn isotopes calcu- 
lated by the Skyrme-HFB theory with the Green’s function 
method (filled circles and solid lines), in comparison with the 
results by the discretized method (open triangles and dashed 
lines). 


clear mass have attracted great interests in the past years 
and lots of efforts have been done, such as the works in 
Refs. [106-112]. 


In the exotic nuclei, diffuse density distributions in the 
coordinate space is often observed. Thus, in Fig. 6, to ex- 
plore the exotic structures in the (a) Ca, (b) Ni, (c) Zr, 
and (d) Sn isotopes, we also plot the neutron density 
distributions p(r), where the solid lines are those ob- 
tained by the Green’s function method in comparison 
with the results by the box-discretized method denoted 
by the dashed lines. It is a global trend that the neu- 
tron density distributions are extended further with the 
increasing neutron number. The shell structures influ- 
ence the density distribution significantly, i.e., compared 
with the bound nuclei, the density distributions of the 
neutron-rich nuclei in the (a) Ca, (b) Ni, (c) Zr, and 
(d) Sn isotopes with the neutron number exceeding the 
neutron closure N = 28,50, 50,82 are much more extend- 
ed, which is consistent with the behaviors of the rms radii 
plotted in Fig. 5. In addition, compared with the Ca and 
Zr chains, the Ni and Sn chains exhibit more diffuse dis- 
tributions in densities, which can be explained by their 
small two-neutron separation energies Son in a large mass 
range plotted in Fig. 1. In the case of the Zr isotopes, the 
density distributions are relatively localized. In combina- 
tion with the large So, plotted in Fig. 1, we are inclined 
to believe the absence of halos in Zr isotopes. Howev- 
er, taking the RCHB theory with the NLSH parameter 
set [113] and the continuum Skyrme-HFB theory with the 
SKI4 parameter set [74, 114], giant halos in Zr isotopes 
have been predicated. In all isotopes, compared with the 
box-discretized method predicting nonphysical sharp de- 
creases of densities at the space boundary, Green’s func- 
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FIG. 6: (Color online) Neutron density p(r) for the (a) Ca, 
(b) Ni, (c) Zr, and (d) Sn isotopes calculated by the 
Skyrme-HFB theory with the Green’s function method (sol- 
id lines), in comparison with the results by the discretized 
method (dashed lines). 
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FIG. 7: (Color online) Contributions of different partial waves 
to the total neutron density, pi;(r)/p(r), as a function of 
the radial coordinate r for the neutron-rich nuclei (a) °Ca, 
(b) ®°Ni, (c) '°°Zr, and (d) '“Sn. The shallow regions are for 
the nuclear density distributions p(r), which are rescaled by 
multiplying a factor of 5. 


tion method can describe well the extended density dis- 
tributions, especially for the very neutron-rich isotopes. 
Besides, the densities obtained by the Green’s function 
method can be independent on the space sizes as dis- 
cussed in Ref. [77], which is determined basically by the 
proper boundary conditions of the bound states, weak- 
ly bound sates and the continuum employed when con- 
structing the Green’s functions. 

To explore the contributions of different partial waves 


on the extended density distributions in Fig. 6, tak- 
ing the neutron-rich (a) ®*Ca, (b) °®Ni, (c) 1?°Zr, and 
(d) !“4Sn as examples, we plot in Fig. 7 the composition- 
s pij(r)/p(r) as functions of the radial coordinate r. It 
can be clearly seen that outside the nuclear surface refer- 
ring to the right boundary of the shallow regions of the 
total nuclear density distributions, it’s the orbitals locat- 
ed around the Fermi surface that play the main role of 
the density distributions. For example, in the neutron- 
rich ©4Ca, the partial waves pi/2, f5/2, 99/2, 81/2, and 
ds/2 contribute a lot for the total density in the area 
of 5 fm < r < 15 fm. Those levels are located within 
~ 5 MeV around the Fermi surface as shown in Fig. 8, 
where the neutron canonical single-particle levels as well 
as the occupation probabilities are presented. When go- 
ing further in the coordinate space with r > 15 fm, the 
contributions of the partial waves s,/2 and d5/2 with s- 
mall angular momenta increase evidently while others re- 
duce. In the case of *°Ni, the single-particle levels 2d5 2, 
381/2, and 2d3/2 are located above the neutron shell of 
N = 50 and close to the Fermi surface, playing the main 
role of the neutron density distribution. Although the 
positive state 1g7/2 is also very close to the Fermi sur- 
face, it does little contributions for the density in large 
coordinate space due to the large centrifugal barrier. For 
the nucleus !2°Zr with the neutron number very close to 
the closure of N = 82, the single-particle levels between 
the closures N = 50 and N = 82 including 2d5/2, 197/2, 
381/2, 2d3/2, and 1h11/2, contribute a lot for the densities 
in large coordinate space. The insignificant occupation of 
the positive state 2f7/2 leads to little contribution to the 
density. Regarding to the neutron-rich '“4Sn with the 
neutron number exceeding the closure of N = 82, the 
weakly bound single-particle levels 2f7/2, 2fs/2, 3p3/2; 
and 3p;/2 with small angular momentum play the key 
role for the extended density distributions in the large 
coordinate space. According to those analysis, we can 
conclude that the single-particle levels around the Fer- 
mi surface and particularly the waves with low angular 
momenta are the main cause for the extended nuclear 
density. 

In Fig. 8, the particle occupation probabilities v? on 
different canonical levels eX" are presented and denoted 
by the length of the lines. Without pairing, the values 
of the occupation probabilities v? should be either one 
or zero separated by the Fermi surface. With the effect 
of pairing, the nucleons occupying the levels below the 
Fermi surface can be scattered to higher levels, which 
results in the occupations of the weakly bound states 
above the Fermi surface and even levels in the continuum. 
In the neutron-rich nuclei ®4Ca, ®°Ni, 12°Zr, and 174Sn, 


the number of neutrons N) = 5 (2j +1)uz scattered 
Ek> An 
above the Fermi surface Ay are 4.33, 2.51, 0.159, and 
0.173, respectively. In the case of ®4Ca, the weakly bound 
single-particle 1gg/z contribute around 3.7 neutrons. 
To explore the effects of pairing, we plot in Fig. 9 the 
number of neutrons N) scattered above the Fermi surface 
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FIG. 8: (Color online) Neutron canonical single-particle lev- 
els in the (a) ĉtCa, (b) 8°Ni, (c) 1°Zr, and (d) "4Sn calcu- 
lated by the Skyrme-HFB theory with the Green’s function 
method. The blue dashed lines represent the Fermi surface. 
The occupation probabilities v? for the canonical orbitals are 
proportional to the lengths of the lines (in logarithmic scale). 


for the (a) Ca, (b) Ni, (c) Zr, and (d) Sn chains obtained 
by the Skyrme-HFB theory with the Green’s function 
method. It can be seen clearly that substantial neutrons 
have been scattered from the single-particle levels be- 
low the Fermi surface to the weakly bound states above 
the Fermi surface and even levels in the continuum due 
to the pairing, especially in the nuclei with the neutron 
number filling the half-full shells. Besides, very obvious 
shell structure could be observed. With the number of 
neutrons reaching a magic number, i.e., N = 28,40 in 
Ca isotopes, N = 40,50 in Ni isotopes, N = 50,82 in Zr 
isotopes, and N = 82,126 in Sn isotopes, N) is almost 
zero due to the absence of the pairing for the closed-shell 
nuclei. Besides, at the points of N = 32 in the Ca iso- 
topes, N = 54,68 in the Zr isotopes, and N = 88 in the 
Sn isotopes, very small numbers of neutron N) are ob- 
tained, which indicate weak pairing in those nuclei. On 
the other hand, it can also predict the possible existence 
of sub-shells and new magic numbers. Furthermore, the 
evolution of N is basically consistent with the trend of 
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FIG. 9: (Color online) Number of the neutrons Ny occupying 
the single-particle levels above the Fermi surface An as a func- 
tion of the mass number A for the (a) Ca, (b) Ni, (c) Zr, and 
(d) Sn isotopes calculated by the Skyrme-HFB theory with 
the Green’s function method. 
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FIG. 10: (Color online) Number of neutrons No occupying in 
the continuum (above the threshold £ = 0 MeV) as a function 
of the mass number A for the (a) Ca, (b) Ni, (c) Zr, and 
(d) Sn isotopes calculated by the Skyrme-HFB theory with 
the Green’s function method. 


pair energy —Fpair in Fig. 3. 

In Fig. 10, we further investigate the number of neu- 
trons occupying the continuum with the single-particle 
energies € > 0 MeV, i.e., No = 5 (2j +1)v?. Com- 

Ek>0 
pared with the number of neutrons N) occupying the 
levels above the Fermi surface, the number of neutrons 
occupying the continuum No is reduced drastically. For 
example, in the Ca chain, No is less than one in all iso- 
topes except °?Ca. In Ca, the single-particle level 1gg/2 
displays as a low-lying canonical positive state in the con- 
tinuum with energy of € = 0.198 MeV and a high occupa- 


tion probability of v? = 0.197, resulting in almost 1.979 
neutrons occupying it. But in the neighboring ©°Ca, very 
small occupation probability v? = 0.015 of 1gg/2 is ob- 
tained while in ®'Ca, the 19/2 state drops to be a weakly 
bound level with the energy of e = —0.072 MeV. After 
removing the neutron contribution from 1gg/2 in °Ca, 
only 0.38 neutrons are in the continuum, which is denot- 
ed by an empty circle in panel (a). Except the Sn chain, 
the shape of No is very close to those of pairing energy 
and pairing gap. Although the case of Sn chain becomes 
very complex, we can still observe the shell structure at 
N = 82 and N = 126 where Np is almost zero. 


V. SUMMARY 


In this work, the exotic nuclear properties of neutron- 
rich Ca, Ni, Zr, and Sn isotopes are examined systemat- 
ically taking the continuum Skyrme-HFB theory in the 
coordinate space formulated with the Green’s function 
method, in which the pairing correlations, the couplings 
to continuum, and the blocking effects for the unpaired 
nucleon in odd-A nuclei are treated properly. 

Firstly, both the two-neutron separation energies Soy 
and one-neutron separation energies Sp are calculated, 
which are consistent with the experimental data. Signifi- 
cant differences exist for the drip lines determined by Soy 
and Sn. In the Ca, Ni, Zr, and Sn isotopes, the drip line 
nuclei are ©’Ca, ®°Ni, 1?3Zr, and !”’Sn judging by Son 
while ©°Ca, °°Ni, 1??Zr, and !48Sn according to Sn. Due 
to the absent contribution of pairing energy by the sin- 
gle unpaired odd neutron, the neutron pairing energies 
—Epair of the odd-A nuclei are around 2 MeV smaller 
in comparison with those of the neighboring even-even 
nuclei. This also explains why the drip lines determined 
by Sa» are much shorter than those by Sən. In addition, 
from the fluctuation trends of the pairing energy, the 
traditional neutron magic numbers are clearly displayed, 
i.e., N = 28,40 in Ca isotopes, N = 40,50 in Ni iso- 
topes, N = 50,82 in Zr isotopes, and N = 82,126 in Sn 
isotopes. 

Secondly, to explore the possible halo structures in the 
neutron-rich Ca, Ni, Zr, and Sn isotopes, the neutron 
single-particle structures, the root-mean-square radii, 
and the density distributions are investigated. In the 
neutron-rich Ca, Ni, Sn nuclei, especially the weakly 
bound nuclei close to neutron drip line, the rms radii 
have sharp increases with significant deviations from the 
traditional r « A'/® rule. Besides, very diffuse spatial 
density distributions are also observed in those nucle- 
i, which reflects the possible halo phenomenon therein. 
By analyzing the contribution of different partial waves 
to the total density, we found that the orbitals located 
around the Fermi surface and particularly those with low 
angular momenta are the main cause for the extended 
nuclear density and large rms radii. 

Finally, the number of halo nucleons which could re- 
flects the effects of pairing are discussed. Two different 


numbers of neutrons are defined, i.e., N) occupying the 
single-particle levels above the Fermi surface An and No 
occupying the continuum. We found that the evolutions 
of N) and No with the mass number A are basically con- 
sistent with the trend of pairing energy —Epair, which 
supports the key role of the pairing correlations in the 
halo phenomena. 
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